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We employ the concept of a dynamical, activity order parameter to study the Ising model in 
a transverse magnetic field coupled to a Markovian bath. For a certain range of values of the 
spin-spin coupling, magnetic field and dissipation rate, we identify a first order dynamical phase 
transition between active and inactive dynamical phases. We demonstrate that dynamical phase- 
coexistence becomes manifest in an intermittent behavior of the bath quanta emission. Moreover, 
we establish the connection between the dynamical order parameter that quantifies the activity, and 
the longitudinal magnetization that serves as static order parameter. The system that we consider 
can be implemented in current experiments with Rydberg atoms and trapped ions. 
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I. INTRODUCTION 

The remarkable progress in the control of ultra cold 
atomic gases and trapped ions has opened a new door for 
studying dissipative many-body quantum systems. Very 
recently it was shown that a carefully designed dissipative 
dynamics arising from an engineered heat bath can lead 
to the formation of pure and coherent many-body quan- 
tum states [IH5]. Moreover, it was demonstrated that 
open many-body systems with competing dissipative and 
coherent interactions possess a rich phase structure |BJ- 
I12j . However, despite their intrinsic dynamical nature, 
these phases are often classified by means of equilibrium 
order parameters, such as particle densities or spatial cor- 
relation functions. 

In this work we pursue a complementary route which 
aims at describing dynamical phases in terms of strictly 
dynamical order parameters. The approach that we will 
use is called thermodynamics of trajectories. It has al- 
ready proven to be useful for the study of classical many- 
body systems displaying complex cooperative dynamics 
such as glasses JT3J HI] . Recent theoretical work [T3] 
has adapted this approach for simple open quantum sys- 
tems and has shown that dynamical phase behavior can 
be uncovered by means of an activity order parameter 
[13l [TBI [17] , i.e. an observable that counts the emission 
of quanta from an open quantum system into its environ- 
ment (events often referred to as quantum jumps). 

Here, we show that the thermodynamics of trajecto- 
ries approach can also be used to gain insights into the 
dynamical behavior of interacting many-body quantum 
systems. We illustrate this by studying a quantum Ising 
model in a transverse magnetic field subject to Marko- 
vian dissipation that couples to individual spins (cf. Fig. 
[l}. Beyond the fact that Ising models like the one we 
study here serve as paradigmatic examples of many-body 
systems, the system at hand can be implemented with re- 
cently developed techniques in experiments with trapped 
ions US] or Rydberg atoms [H [H [H] ■ 

We show that depending on the experimental param- 
eters the system is either found in a specific dynamical 
phase or at coexistence conditions of two such phases. 
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FIG. 1. (Color online) (a) Schematic of the open Ising model 
in a transverse field. Internal states of trapped atoms or ions 
are described by spin-1/2 degrees of freedom on a lattice. 
The interaction between nearest neighbors has strength V, 
and there is an applied transverse field f2 (e.g. the Rabi fre- 
quency of a laser). Interaction with the radiation field leads 
to the incoherent emission of bath quanta (e.g. photons) with 
rate k, which are detected and counted, (b) Quantum jump 
trajectory in a one-dimensional version of (a). Emitted pho- 
tons are temporally (and spatially) resolved, so each point 
indicates where and when a quantum jump event took place. 
This particular trajectory shows intermittency, which mani- 
fests the coexistence of an active and an inactive dynamical 
phase (see text). 



We demonstrate that this phase coexistence is accompa- 
nied by strong fluctuations in the activity that become 
manifest in pronounced intermittency in the emission of 
the bath quanta. The thermodynamics of trajectories 
approach also sheds light on the intermittency that has 
recently been theoretically found in an open fully con- 
nected spin model [12] . 

we 



The paper is organized as follows. In Section II A 



outline the thermodynamics of trajectories approach for 
quantum systems. In particular, we establish the rela- 
tion between the dynamical phases, as classified by their 
activity, and equilibrium phases characterized by static 
observables (Sec. II B). In addition, we show that such a 
direct connection does, in general, not hold between dy- 
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namical and static fluctuations. Following these general 
considerations, we study the dynamical behavior of the 
dissipative quantum Ising model in a transverse field us- 
ing the thermodynamics of trajectories approach in Sec- 
tion III After briefly discussing a possible realization of 



this model with highly excited atoms (Sec. Ill A I we in- 
vestigate the dynamics of the system on the level of mean- 
field theory (Sec. IIIB) as well as by using the numer- 

and 
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ical techniques of exact diagonalization (Sec. 
Quantum Jump Monte-Carlo simulations (Sec. 
focus, in particular, on the regime in which the model 
displays strongly intermittent behavior in the emission 
of bath quanta. We show that this is a consequence of 
the system being at a first order coexistence point between 
two dynamical phases. 



II. THERMODYNAMICS OF TRAJECTORIES 

A. Dynamical order parameter and ensembles of 
quantum trajectories 

Let us consider a quantum many-body system com- 
posed of N particles coupled to a Markovian bath. The 
state of the system is described by the density matrix p 
and its time-evolution governed by the Master equation 
dtp = W(p) with super-operator, 

N N 

w(.) = -i [h, .] + J2 Ju • 4 - - E { J Uu,*} , (i) 



and system Hamiltonian H [21]. The coupling to the 
thermal bath is assumed to take place at the single par- 
ticle level, and is mediated by the Lindblad operators 
J„ [25 . These operators induce sudden changes of the 
state of the system and, therefore, describe spontaneous 
emission events of quanta into the bath. 

We are interested in the statistics of these emission 
events and their connection to the quantum dynamics of 
the system. To illuminate this we employ the thermody- 
namics of trajectories approach, which will allow us to 
gain immediate insights into the emission characteristics 
of bath quanta. The idea, which for quantum systems 
is outlined in Ref. [15] . is to consider time records of 
emission events — so-called quantum jump trajectories — 
and to study ensembles of trajectories in a manner anal- 
ogous to how one studies ensembles of microstates (or 
configurations) in equilibrium statistical mechanics. A 
trajectory of length (i.e. observation time) t contains K 
emission events; see e.g. Fig. [ljb). The probability for 
such a trajectory to occur is given by px (t) = Tr p( K 1 (t) 
where (<) is the projection of the density matrix p 
onto the subspace in which K emission events have taken 
place [22] ■ The generating function of the probability dis- 



tribution px (t) is defined as 

Z=jre-s K PK (t)=jre-° K TrpW(t) 

oo 

^Tr^e- s V^(<)^Trp s (i). (2) 



K=0 



Here p s (t) is the Laplace transformed density matrix 
which evolves according to the generalized Master equa- 
tion d tPs = W s ( Ps ) = W(p.) + V s ( Ps ) MM where 



V s (.) = (e- s - 1 



This generalized Master equation is not probability con- 
serving, and corresponds to non-physical dynamics where 
the probability of quantum jump trajectories with K 
events are biased by a factor e~ sK . This means that for 
negative (positive) s trajectories with more (less) emis- 
sion events than the average are more likely to occur. The 
physical (probability conserving) dynamics takes place at 
s = 0, where the generalized Master operator coincides 
with ([I]). However, as we show below, information about 
the behavior of the system in the vicinity of this physical 
point (i.e. at s ^ 0) can be crucial for the understanding 
of its emission dynamics [HI [25] . 

We are interested in dynamical properties in the sta- 
tionary regime. Therefore, we consider observation times 
long enough for all initial transient effects to have become 
negligible. For these long times the generating function Z 
acquires a large-deviation form [26, 27 . Using the spec- 
tral decomposition of the generalized Master operator we 
can write in the limit of long times Z = Tr p s (t) — > e' e W 
where 8(s) is the eigenvalue of W s with the largest real 
part. The crucial idea behind the thermodynamics of tra- 
jectories approach is to interpret Z as a partition function 
and of 9(s) as a free-energy, i.e. to ascribe a real physi- 
cal meaning to these quantities. This allows us to view 
the (s-dependent) mean emission rate (k)(s) = (K)(s)/t 
of bath quanta as a dynamical order parameter which 
we call the activity. The activity can be written as the 
derivative of the free-energy with respect to the conju- 
gate field s: (k)(s) = —d s 8(s). Likewise, higher mo- 
ments of the emission statistics are encoded in higher 
derivatives of 6(s). Of particular interest here is Man- 
del's Q-parameter, Q(s) — ~d^8(s) /d s 9(s) — 1, which 
quantifies the deviation of the emission statistics from a 
Poissonian distribution. 

We would like to remark that technically the ther- 
modynamics of trajectories approach bears similarities 
to that of "Full Counting Statistics" (FCS) [23 ESI- 
However, despite the similarities to FCS we approach 
the counting problem from a somewhat different angle. 
We regard the counting field s as the conjugate field 
to a dynamical order parameter (the event count K). 
This very perspective allows to construct a theoretical 
framework similar to equilibrium statistical mechanics 
for the analysis of ensembles of stochastic trajectories. 
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For example, just like in equilibrium statistical mechan- 
ics, non- analytical points of the dynamical free-energy 
9(s) will determine transitions between dynamical phases 
[T51 1251 [5U] . Of particular importance are of course non- 
analyticities that occur at or very near the physical point 
s = 0, as they strongly influence the dynamical behavior 
of the system. 

Before we further illustrate this point by investigating 
the emission statistics of a particular model system, let 
us use the approach outlined above to assess possible 
connections between dynamical and static quantities. In 
particular, we will show that the activity can directly be 
related to the expectation value of a static observable. 
We will also see that, in contrast, the same is in general 
not true for fluctuations (and higher moments) of the 
dynamical order parameter. 



This is in fact known from strongly interacting classi- 
cal systems. There, dynamical fluctuations are not nec- 
essarily simple manifestations of static fluctuations. A 
notable example is glasses: there, the thermodynamics 
can essentially be trivial, while the corresponding dy- 
namics can be very complex. For example, Refs. [321 [14] 
show that both idealized lattice models and realistic liq- 
uid models, when explored by means of techniques sim- 
ilar to those employed here, display phase transitions in 
dynamical trajectories which have no static counterpart. 
To uncover these dynamical transitions it is necessary 
to go beyond purely static treatments and consider the 
statistical properties of trajectories. 

In the following we will apply the general ideas of the 
thermodynamics of trajectories approach to study the 
dynamical behavior of a specific open many-body quan- 
tum system. 



B. Connection between dynamic and static 
quantities 

As a consequence of the large deviation principle, the 
dynamical free-energy for long times can be written as 
0(s) = dthxZ = d t lnTrp s = (Tr d t p s )/(Tr p.). Using the 
equation of motion for p s , we get Tr9tp s = TrV s (p s ), 
where we have exploited that TrW(p s ) = 0. Combining 
these results yields 



9(a) = (e~ 



-l)1t^4J v R(a)\ 



(3) 



where we have introduced the normalized density matrix 
R(s) — ps/Trps, which becomes stationary in the long 
time limit. Furthermore R(s) coincides with the station- 
ary density matrix of the system at the physical point, 
R(s = 0) = p. Differentiation of Eq. ^ with respect to s 
reveals the relation between the dynamical order param- 
eters and the expectation value J^J V ) at s = 0, 



(k) (0) = -d s 6(s)\ s=0 = TrJ2jU 



p. 



(4) 



A further differentiation of Eq. |3| shows that this di- 
rect connection between dynamical and static quantities 
does not hold for higher moments. In particular, tempo- 
ral fluctuation of emission records (quantum jump tra- 
jectories) are in general not solely determined by static 
spatial correlations. Assessing higher moments of the 
event count distribution, therefore, requires knowledge 
of the dynamical phases of the system in the vicinity of 
the physical point, i.e. for s ^ 0. This can, e.g. been 
seen when evaluating the Q-P arame ter, 



Q(Q) = _ 2 ^E,iJW(4. 



(5) 



Calculating Q requires the evaluation of the first deriva- 
tive of R(s) with respect to the field s, i.e. information 
about R(s) for s ^ 0. 



III. DISSIPATIVE ISING MODEL IN A 
TRANSVERSE MAGNETIC FIELD 

We consider a system of N spin-1/2 particles arranged 
on the sites of a regular lattice with lattice spacing a and 
coupled to a thermal bath with zero temperature. The 
Hamiltonian H is that of a quantum Ising model in a 
transverse field [3T] 



(6) 



with Q characterizing the transverse field strength and 
V being the interaction strength between adjacent spins 
((//, v) indicates that the sum is only over nearest neigh- 
boring pairs). The Lindblad operators are given by 



Ju = \/kS^ = y/K[Sx"' — iSy"'] with decay rate K. This 
model is sketched in Fig. [lja) . 

A natural static order parameter for this system is the 
magnetization m = N~ x ^2 V (S^), which has also been 
used to characterize dynamical properties in Refs. 
ITI2] . Using Eq. Q we immediately see that the static 
order parameter is indeed proportional to the dynamical 
order parameter, 



Mi 



(k) (0) = kN 



(7) 



The static magnetic susceptibility, however, is not di- 
rectly related to quadratic fluctuations in the activity, as 
the former only captures static spatial correlations, while 
the latter also includes correlations in time. 



A. Realization with Rydberg atoms 

Experimentally the dissipative Ising model can - to a 
good degree of approximation - be achieved with elec- 
tronically excited ultracold alkali atoms confined to a 
lattice. Here, the state \\.) is identified with the atomic 
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FIG. 2. (Color online) Mean-field phase diagram. Colored 
regions have a unique steady state solution. In the grey re- 
gion, delimited by the two spinodal lines vi and V2 (see text), 
two steady state solutions with different magnetization exist. 
These solutions merge at (Q*;V*) = (2k; (3\/3/z)k). The 
small panels show the behavior of the free-energy 8(s) (solid) 
and its first derivative (k)(s) (dotted) in the vicinity of s = 0. 
While in the region with unique steady state both functions 
are smooth, the dynamical order parameter (k) (s) has a jump 
at s = for parameter values taken in the grey region. This 
indicates the existence of a first order dynamical phase tran- 
sition. 



ground state and the state |f) is an electronically ex- 
cited Rydberg nS-state. When two atoms on the v- 
th and ^-th lattice site (with position vectors r v and 
r^, respectively) are excited simultaneously they inter- 
act via a van-der-Waals potential of the form V„ p — 
\ Tv C -%^ |t)„(t| ® 11% (tl- Here C 6 is the dispersion coef- 
ficient which characterizes the interaction strength. Due 
to the quick decay of the interaction as a function of 
the distance one can replace the van-der-Waals potential 
by a nearest-neighbor interaction [35] ■ When Rydberg 
states are excited by a laser of Rabi-frequency f2 and 
detuning A with respect to the energy difference of the 
transition \l) O |t), the Hamiltonian for the f-th atom 
reads h v — QSx + AjS* . Combining this with the in- 
teraction one finds that by choosing A = — V = —C^/a 6 
one arrives at Hamiltonian Q . The dissipative dynamics 
is realized naturally as Rydberg states decay radiatively. 
For nS'-states of alkali metal atoms this decay takes place 
predominantly to the lowest P-state (on the time-scale 
of a few microseconds) and subsequently to the ground 
state. The last decay is very fast (tens of nano-seconds) 
such that it can be considered instantaneous. Within this 
approximation we find the dynamics of the system to be 
governed by the Master operator 0. The activity can 
then be directly monitored by detecting photons that are 
emitted during the decay from the lowest P-state to the 
ground state. 



B. Mean- field approximation to statics and 
dynamics 

The static properties of our open Ising model, includ- 
ing the nature of its structural phases and the transitions 
between them, are given by its stationary density matrix. 
This is obtained from the eigenstate(s) of the operator W 
with zero eigenvalue [21]. Finding the exact eigenstate(s) 
is a difficult task, but we can learn much from a mean- 
field approximation [TTJ [T^] . 

Since the dynamical order parameter is directly pro- 
portional to the static one [cf. Eq. Q], determining the 
activity is equivalent to calculating the magnetization. 
In mean-field approximation this can be done in a stan- 
dard manner by approximating the many-body density 
matrix p as a product of single-particle density matrices, 

P ~ ®f=i and solving Tr 2 ...iv W (®f =1 = 

[3"3"] . The partial traces are taken over the degrees of free- 
dom of all spins but one. Parameterizing the single parti- 
cle matrices as r<"> = t/2+aS^ +/3S^ ) +mSi" } this 
procedure leads to the following self-consistency equation 
for the magnetization m: 




= 8z 2 V 2 m 3 + 4z 2 V 2 m 2 + 2(k 2 + 2n 2 )m + k 2 , (8) 

where z denotes the coordination number of the lattice. 
The analysis of this third order polynomial shows that 
below a critical interaction strength V* = (3\/3/z)k the 
self-consistency equation ^ has only one stable real so- 
lution. Above V* one finds, depending on the value of f2, 
either one or two stable real solutions. The two spinodal 
lines [34] separating the regions of unique and multiple 
real solutions for m are given by 
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valid for V/k, Q/k > 1. 

The corresponding phase diagram is shown in Fig. [2] 
In the colored regions there is a unique steady state. In 
the grey domain two stable solutions exist, one with small 
magnetization hia ~ and one with large negative mag- 
netization mi ~ —1/2. By virtue of Eq. ^ we can now 
conclude that the activity (k) is small (large) in the blue 
(red) regions. To understand the dynamical behavior of 
the Ising system in the grey region we have to explore 
the behavior of 9(s) in the vicinity of s = 0. 

For this, we expand 9(s) given in Eq. Q to first or- 
der in s and find that 9(s) -sY, v TWj J v R(0) = 
—s(k)(0) = —sKN(m+ 1/2). Taking the magnetization 
obtained from the mean-field calculation we see that in 
the colored regions of Fig. [2] we obtain a 6(s) whose first 
derivative is smooth (see top left and bottom right panel) . 
In the grey region one has to compare the values of 9(s) 
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that one obtains for the two solutions wia and mi, choos- 
ing the one that maximizes 0(s), hence 

s > : 9(s) = -skN (mi + - 



s < : 9(s) = -skN yrriA + ^ 

Since mA ^ mi, the slope of 9(s) changes at s = 0, 
causing a non-analyticity of 9(s) and a jump of its first 
derivative, the activity (see top right panel of Fig. [2]). The 
grey area in the phase diagram can thus be regarded as 
a coexistence region of two dynamical phases: an active 
and an inactive one. Trajectories in the active phase 
are dense in quantum jumps and are characterized by a 
large (k); quantum jumps are scarce in trajectories in the 
inactive phase, which is characterized by a small (k). The 
discontinuity of (k) at s = indicates that the transition 
between active and inactive phases is of first order. 

To illuminate the implications of these results on the 
emission dynamics of the system let us consider a thermo- 
dynamic analogy: a fluid system at the transition point 
between a high density liquid phase and a low density 
vapor one. Here, a small change in pressure will cither 
select the liquid or the vapor, which are distinguished 
by their average specific volumes, confirming that the 
system is at a point of phase coexistence. Therefore, 
the ability to vary the pressure yields important infor- 
mation about the nature of the phase transition. As a 
further consequence of the first-order coexistence one ex- 
pects sharp interfaces between the phases. These inter- 
faces are not a property of either phase, but of the fact 
that there exists a "surface tension" between the phases. 

In our dynamical case the field s works in exactly the 
same way as the pressure in our liquid-vapor analogy - 
it selects (depending on its sign) an active (bright) or 
inactive (dark) phase. This allows us to uncover a first 
order dynamical phase transition and again shows the 
necessity of studying the system away from s — 0. Fur- 
thermore, we expect to observe sharp interfaces between 
the dynamical phases. Since quantum jump trajectories 
of many-body systems live in space and time, the inter- 
faces between distinct dynamical phases can be tempo- 
ral. The mean-field results, therefore, indicate an inter- 
mittent emission pattern, for system parameters chosen 
from the grey region of the phase diagram shown in Fig. 
[2] This region of phase coexistence ends at the point 
(ft*; V*) — (2k; (3\/3/z)k), reminiscent of the static crit- 
ical point beyond which liquid-vapor coexistence is pos- 
sible [33. 

The mean-field approximation above disregards the fi- 
nite range of the spin-spin interactions. Such a treat- 
ment would be accurate for a fully connected problem 
(i.e., every spin interacts with all other spins with the 
same strength), as e.g. studied in [T2]. There, one might 
indeed anticipate strongly collective dynamics. However, 
the spin system studied in this work exhibits local inter- 
actions, and, therefore, it is not evident a priori (particu- 
larly in low dimensions) that the collectiveness in photon 
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FIG. 3. (Color online) Dynamical free-energy 9(s) deter- 
mined by exact diagonalization of the generalized Master op- 
erator for the one-dimensional, dissipative Ising model in a 
transverse field for N = 6 spins with periodic boundary con- 
ditions. The plot shows 9(s) for V = 100k at three different 
values of the coupling strength Q corresponding to: an active 
(red) and inactive (blue) dynamical phase and the behavior 
of the dynamical free-energy under conditions of phase coex- 
istence (black line). 



emission or the intermittence predicted by a mean-field 
analysis should persist when going beyond this approxi- 
mation, as we will do in the following. 



C. Numerical analysis: Exact diagonalization 

To check whether the mean-field treatment yields the 
qualitatively correct behavior of 9(s) for our system with 
short range interactions, we have performed an exact nu- 
merical diagonalization of the generalized Master opera- 
tor for chains up to N = 7 spins. The largest eigenvalue 
corresponds to the dynamical free-energy 9(s) (see e.g. 
Refs. [HIES]). 

The results for 9(s) for N = 6 spins are shown in Fig. 
|3] as a function of s for V — 100k and different values 
of SI. At Q = 42k (red curve) we find the system in 
the active phase which is indicated by a large slope at 
s = 0. Conversely, at Q = 2k (blue curve) the activity 
is small as the slope of 9(s) at s = is almost zero. In 
the region in between [fl = 28k (black curve)] we find a 
crossover compatible with a smoothed first order phase 
transition at s = 0. This result corroborates the expecta- 
tions for the large N limit obtained from the mean-field 
calculation. In fact the curves shown here are a smoothed 
version of the ones shown in Fig. [2j 

To learn more about the dynamical behavior of the 
system at the physical point (s = 0) we have cal- 
culated the activity (k) = — d s 9(s)\ s= o (lines in left 
graph of Fig. El and Mandel's Q-parameter, Q — 
-d^9(s)\ s=0 /d s 9(s)\ s=Q - 1, (lines in right graph of Fig. 
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FIG. 4. (Color online) Numerically determined activity {k)/n (left) and Mandel's Q-parameter (right) as a function of 
Q/k for the one-dimensional, dissipative Ising model in a transverse magnetic field for AT = 5 ... 7 sites and periodic boundary 
conditions. The lines show the results of an exact numerical diagonalization of the generalized Master operator and the symbols 
the results of Quantum Jump Monte-Carlo simulations (ensembles of 1000 trajectories of length t = 200k -1 obtained by slicing 
a much longer one of length t max — 2 x 10 5 k _1 ). 



[4]) as a function of the transverse field strength for system 
sizes of N = 5 . . . 7 spins. In addition to the results of the 
exact diagonalization approach, Fig. [4] also shows data 
extracted from Quantum Jump Monte-Carlo simulations 
(symbols) that will be discussed in the next subsection. 

The activity - being proportional to the magnetization 
m - shows the expected behavior as function of the trans- 
verse field coupling strength fi. For small O (m k — 1/2) 
the system is in an inactive phase with (k) ss 0. In con- 
trast, for large f2 (m w 0) the system is in an active phase 
and the activity saturates at (k)/n — N/2 as expected 
from Eq. Note that due to the finite size of the sys- 
tem the crossover between these phases is smooth. The 
Q-parameter is close to zero in the inactive (small fl/n) 
as well as in active phase (large Q/k) indicating a Poisso- 
nian distribution of the event counts in both dynamical 
phases. It becomes large and positive in the crossover 
region, the height of the peak growing with the system 
size and its position slightly shifting towards larger val- 
ues of Q/k. This region is where we anticipate dynamical 
phase-coexistence in the thermodynamic limit (t oo, 
(A') — > oo while (K)/t = const < oo). This is an indica- 
tion (but no proof) of the mixing of the two dynamical 
phases as would be expected in a region of first-order 
phase coexistence. 

To obtain further evidence that the phase mixing sce- 
nario does indeed hold, one needs further information 
about the distribution px of the event counts. This can, 
in principle, be done by determining higher moments of 
Pk using higher derivatives of 0(s). This, however, will 
quickly get intractable in practice, so we will follow an 
alternative route to show that the picture of the mix- 
ing of two dynamical phases does describe the emission 
dynamics of the system. 



D. Numerical analysis: Quantum Jump 
Monte-Carlo simulations 

To gain a direct insight into the emission characteris- 
tics, we numerically study an unravelling of the full dy- 
namics under the Master operator ([I]) for spin chains of 
up to N = 12 sites using Quantum Jump Monte Carlo 
(QJMC) simulations [35]. To connect with potential ex- 
periments we can interpret a single quantum jump tra- 
jectory of the numerical simulation as a time series of 
emission events recorded by a detector with temporal 
and (for illustration purposes here) also spatial resolu- 
tion, see Fig. [IJx For each set of parameters (Q, V, k) 
this stochastic dynamics generates an ensemble of quan- 
tum jump trajectories. 

In Fig. [4] we compare the mean activity and the Q- 
parameter obtained from QJMC simulations for system 
sizes of N = 5 ... 7 with the results from the exact di- 
agonalization discussed above and find good agreement 
between them. However, while the activity can be ob- 
tained very accurately, the error bars for Q, particularly 
at intermediate values of Q, are relatively large. This 
signifies the strongly fluctuating character of the emis- 
sion dynamics in the crossover region and shows that in 
order to capture the dynamical behavior of the system 
in a possible experiment very long trajectories will be 
required. Faithfully determining higher moments of the 
event count distribution will quickly become intractable. 
However, further information about the counting statis- 
tics can still be obtained by resorting to an alternative 
way to analyze the emission records. To demonstrate 
this, we have simulated the emission statistics of a chain 
of N = 12 spins and analyzed the data by determining 
the distribution of activities and studying the distribu- 
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FIG. 5. (Color online) Analysis of data from Quantum Jump Monte Carlo simulations of a one-dimensional open Ising model in 
a transverse field with N = 12 spins (k = 0.1). (a) Activity distributions of trajectories of length t = 10k -1 (green) and 50k _1 
(transparent with solid lines), for three different values of fi/re for fixed V = 100k. These ensembles of trajectories are obtained 
by slicing longer trajectories of length t max = 10 4 k -1 . We also show sample trajectories at the corresponding conditions, in 
which the emission events are site resolved, (b) Waiting time distribution p(r) for Q = 25k and V = 100k (crosses in green). 
The lines represent exponential waiting time distributions with time constants 1/&A (solid red line) and 1/fci (dashed blue line), 
(c) Q T parameter as a function of the transverse magnetic field for V = 100k extracted from the waiting time distribution. The 
insets show samples of typical trajectories in each parameter regime. The maximum of Q T coincides with highly intermittent 
dynamics. 



tion of the waiting time between quantum jumps (Fig. 
§. 

The mean-field analysis indicates that for certain com- 
binations of (f2, V, k) we should encounter a dynamical 
phase transition of first order. The hallmarks of such a 
transition would be a bimodal distribution of the activity 
order parameter and sharp interfaces at coexistence con- 
ditions. We generate ensembles of trajectories of length 
t by slicing much longer trajectories of length i max ^> t. 
The time t can vary, but has to be long enough so that 
the dynamics within each phase are captured, but not too 
long compared to the persistence or survival time within 
each phase which is finite due to the finite size of the sys- 
tem. In Fig. [5ja) we show activity histograms at several 
values of O/k for fixed V and two different values of the 
time slicing t = 10k -1 and 50k -1 (note that the results 
for both cases are consistent). For small (large) values 
of f2 the activity is mostly on the inactive (active) side, 
as expected from the mean-field analysis. These distri- 
butions show bimodality at intermediate values of Q/k: 
one peak corresponds to active and the other one to in- 
active dynamics. When this occurs, the corresponding 
trajectories are highly intermittent, showing prolonged 



periods of activity and prolonged periods of inactivity, 
delimited by sharp temporal interfaces as expected from 
a first order transition scenario. 

Further evidence of this behavior is provided by the 
waiting time distribution p(r) of the emission events 
[see Fig. [5](b)] which shows the existence of two distinct 
timescales corresponding to the typical waiting times 
within each phase that govern the dynamics of the sys- 
tem. These two timescales are 1/fci and 1/&A, which 
directly relate to the activities of the coexisting inactive 
and active phases, that are given by k\ and k&, respec- 
tively. 

The results shown in Fig. [5]ja) and (b) suggest that 
at the coexistence region the dynamical fluctuations of 
the open Ising model are dominated by the switching 
between active and inactive behavior, rather than fluc- 
tuations within each of these phases. Disregarding intra- 
phase fluctuations completely we can approximate the 
waiting time distribution by the sum of two exponentials, 

, , ak?e~ klT + (1 - a)k\e- kAT 

p ^ = T~Tn \& ' 

ak\ + (1 — a)kA 

i.e. we assume that the active(inactive) dynamics is 



Poissonian with average number of events per unit time 
kA{ki). The lines in Fig. |5jb) represent these two con- 
tributions, where k\ and are extracted from the nu- 
merical simulation. The mixing of the two dynamical 
phases is described by the parameter a — ... 1, which 
determines the point, at which the two lines in Fig.[5]jb) 
cross (in the particular case shown a ~ 0.5). The dy- 
namical behavior in the intermediate region is strongly 
mixed and highly non-Poissonian due to the coexistence 
between two very distinct (close to) Poissonian phases. 

Intermittency in the dynamical trajectories corre- 
sponds to large fluctuations of the times between quan- 
tum jumps. In order to quantify these fluctuations 
we use the function Q T = (t 2 )/(t) 2 — 2 where 
(t a ) = J °° dr T X p{r) stands for the A-th moment of 
p(r). This function reaches its minimum Q T = — 1 if the 
waiting time distribution gives rise to a completely regu- 
lar distribution of the quantum jumps without statistical 
fluctuations. It assumes the value Q T — when p(r) is 
an exponential distribution, that is, if the distribution 
of quantum jumps is Poissonian. When the fluctuations 
of times between jumps are large, it assumes a positive 
value Q T > 0. The behavior of Q T as a function of Q/k is 
shown in Fig. ]5jc) . It does indeed peak for values of the 
parameters for which the activity is bimodal and the tra- 
jectories intermittent. Far from the coexistence region 
Q T drops to values close to zero, suggesting that fluc- 



tuations within each phase are well approximated by a 
Poisson process, at least for the long timescales relevant 
for the discussion here. 



IV. CONCLUSIONS 

We have shown that the concept of a dynamical activ- 
ity order parameter can be successfully applied to char- 
acterize and understand complex dynamical behavior of 
a many-body quantum system. We have established a 
general connection between static observables and the 
dynamical order parameter. For the dissipative Ising 
model studied here we have identified a first order dy- 
namical phase transition, and showed that coexistence of 
two dynamical phases gives rise to pronounced intermit- 
tency of the bath quanta emission. 
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